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Abstract. We examine the climate effects of the emissions 
of near-term climate forcers (NTCFs) from 4 continental re- 
gions (East Asia, Europe, North America and South Asia) 
using results from the Task Force on Hemispheric Transport 
of Air Pollution Source-Receptor global chemical transport 
model simulations. We address 3 aerosol species (sulphate, 
particulate organic matter and black carbon) and 4 ozone 
precursors (methane, reactive nitrogen oxides (NO x ), volatile 
organic compounds and carbon monoxide). We calculate the 
global climate metrics: global warming potentials (GWPs) 
and global temperature change potentials (GTPs). For the 
aerosols these metrics are simply time-dependent scalings of 
the equilibrium radiative forcings. The GTPs decrease more 
rapidly with time than the GWPs. The aerosol forcings and 
hence climate metrics have only a modest dependence on 
emission region. The metrics for ozone precursors include 
the effects on the methane lifetime. The impacts via methane 
are particularly important for the 20 yr GTPs. Emissions of 
NO x and VOCs from South Asia have GWPs and GTPs of 
higher magnitude than from the other Northern Hemisphere 
regions. 

The analysis is further extended by examining the 
temperature-change impacts in 4 latitude bands, and cal- 
culating absolute regional temperature-change potentials 
(ARTPs). The latitudinal pattern of the temperature response 


does not directly follow the pattern of the diagnosed radia- 
tive forcing. We find that temperatures in the Arctic latitudes 
appear to be particularly sensitive to BC emissions from 
South Asia. The northern mid-latitude temperature response 
to northern mid-latitude emissions is approximately twice as 
large as the global average response for aerosol emission, and 
about 20-30% larger than the global average for methane, 
VOC and CO emissions. 


1 Introduction 

The emissions of reactive gases and aerosols can influence 
human and ecosystem health by affecting ozone and par- 
ticulate matter concentrations (HTAP, 2010). They can also 

affect climate through the burdens of ozone, methane and 
aerosols, having both cooling and warming effects. Because 

of the short lifetimes of aerosols (days), ozone (weeks), 
methane (a decade), and their precursors, their climate ef- 

fects are predominantly in the near term (less than 30 yr) 
so we refer to these species as “near-term climate forcers 
(NTCFs)”. They are also often called short-lived climate pol- 
lutants (SLCPs). 

UNEP and WMO (2011) and Shindell et al. (2012) have 
suggested that the mitigation of ozone precursors and black 
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carbon (BC) would be attractive for both air quality and cli- 
mate on a 30-yr timescale, provided it is not at the expense 
of CO2 mitigation. Typical air quality policies target both 
warming and cooling species and tend to have an overall 
detrimental effect on climate (in terms of surface tempera- 
tures). Therefore, it is important to understand how the ef- 
fects of NTCFs vary by location of emissions, when air qual- 
ity policies are considered. This can also be important for 
climate policies that affect both short and long-lived species 
(Berntsen et al., 2006). Many metrics have been proposed to 
compare these effects on climate, but here we consider the 
integrated radiative forcing (RF) using the global warming 
potential metric (GWP) (IPCC, 1990) and the surface tem- 
perature change using the global and regional temperature- 
change metrics (GTP and RTP) (Shine et al., 2005; Shindell 
and Faluvegi, 2010). 

The NTCFs we consider in this paper are sulphate, par- 
ticulate organic matter (POM), black carbon (BC), methane 
and ozone precursors. Some halogenated species have short 
lifetimes and are therefore NTCFs, but we do not consider 
those here. NTCFs with lifetimes longer than the interhemi- 
spheric mixing time (such as methane) are considered rea- 
sonably well mixed, with the concentrations and RF patterns 
independent of the emission location. Species with shorter 
lifetimes such as ozone and aerosols have heterogeneous dis- 
tributions and RF patterns that are dependent on the emis- 
sion location (Fuglestvedt et al., 1999; Berntsen et al., 2005; 
Naik et al., 2005; Fry et al., 2012; Yu et al., 2013). The sur- 
face temperature response does not directly follow the spatial 
details of the RF pattern, rather it smoothes out the pattern 
over scales of '3500 km in the meridional direction and over 
12000 km in the zonal direction (Shindell et al., 2010). 

The RFs and GWPs for the Task Force on Hemispheric 
Transport of Air Pollution (HTAP) Source-Receptor global 
chemical transport model (CTM) simulations (HTAP, 2010) 
were previously documented by Fry et al. (2012) for reac- 
tive gases, and by Yu et al. (2013) for aerosols. In this study 
we have moved further down the chain of climate impacts, 
showing how analytical formulae can be used to relate equi- 
librium RF values to the evolution of global and latitudinal 
temperature changes. 

We use the HTAP RF estimates to calculate GWPs in 
Sect. 3.1. From these, we derive the global mean surface tem- 
perature responses as functions of time for emissions of dif- 
ferent components in different Northern Hemisphere regions, 
and normalise by the corresponding response for CO2 to de- 
rive GTPs in Sect. 3.2. Sections 3.1 and 3.2 follow a similar, 
but not identical methodology to Fuglestvedt et al. (2010). 
One difference being that for the ozone precursor reduction 
cases, we use experiments with fixed global methane con- 
centrations; our calculations are based on RF and changes 
in methane lifetimes rather than the steady-state atmospheric 
concentrations and primary mode time constant that were 
used in Fuglestvedt et al. (2010). Then in Sect. 4, we go 
beyond looking at global mean values for the response and 


use the recently introduced concept of the RTP (Shindell and 
Faluvegi, 2010; Shindell et al., 2012) to quantify the impacts 
of regional pollutant emissions on temperatures in broad lat- 
itude bands. 


2 Methodology 

The basic data needed as input to the GWP and GTP cal- 
culations can be derived from the equilibrium RF responses 
of chemistry transport models to step changes in emis- 
sions. Here we use results from the HTAP study (HTAP, 
2010) in which an ensemble of chemistry transport mod- 
els carried out experiments reducing the emissions of (an- 
thropogenic) ozone or aerosol precursors by 20% in four 
Northern Hemisphere continental regions (North America 
[NA, 15°-55° N, 60° W-125° W], Europe [EU, 25°-65° N, 
10° W-50° E], East Asia [EA, 15°-50° N, 95°-160° E], and 
South Asia [SA, 5°-35° N, 50°-95° E]). Each ozone precur- 
sor was reduced separately in each region. For the aerosol 
experiments, all species were reduced together. The effects 
of regional ozone precursor emission changes (NO x , CO, 
NMVOCs, methane) on ozone and methane concentrations 
were calculated by Fiore et al. (2009). The consequent net 
RF impacts from an ensemble of 1 1 of the CTMs were cal- 
culated by Fry et al. (2012). The changes in RF due to the 
emissions of SO2, primary particulate organic matter (POM) 
and BC were calculated by Yu et al. (2013) from an ensemble 
of 9 models. 

Because of the time needed for methane to equilibrate to 
emissions changes, methane concentrations were fixed glob- 
ally to 1760 ppb (nmol mol -1 ) for all experiments except for 
the methane perturbation experiment in which it was reduced 
by 20 % to 1408 ppb. The impacts of ozone precursor emis- 
sions on global methane and long-term ozone responses were 
calculated analytically from the methane loss by tropospheric 
OH diagnostic reported by each CTM (Fiore et al., 2009; Fry 
et al., 2012). 

Fry et al. (2012) used the NOAA Geophysical Fluid 
Dynamics Laboratory (GFDL) standalone radiative trans- 
fer model (GFDL GAMDT, 2004) to calculate the grid- 
ded, stratospheric-adjusted net radiative fluxes (incoming 
shortwave minus outgoing longwave) at steady state due to 
changes in methane, ozone, and sulphate aerosol. The 3-D 
concentration fields were from the multi-model mean along 
with the mean ± 1 standard deviation. The radiative transfer 
model simulations exclude aerosol-cloud interactions and the 
internal mixing of aerosols. 

Yu et al. (2013) applied normalised direct RFs (NDRFs) 
calculated with the Goddard Chemistry Aerosol Radiation 
and Transport (GOCART) model (Yu et al., 2004) to the 
aerosol optical depths (AODs) from each model. The NDRFs 
were applied on a 2° x 2.5° degree grid for each month of 
the year for each of the three components (SO4, POM, BC). 
No aerosol-cloud interactions were included and the aerosols 
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Table 1. Absolute global warming potentials (AGWPs) for emissions of SO 2 , POM and BC from the 4 HTAP continents, based on Table 6 
of Yu et al. (2013). Units are Wm - ’ yrkg -1 .The AGWPs for aerosols are approximately independent of time horizon. Note these values 
exclude the effects of aerosol cloud interactions and the effect of BC deposition on snow albedo. 



EA 

EU 

NA 

SA 

All 

so 2 

POM 

BC 

—2.9 ± 0.8 x 10 -12 
— 3.7 ± 1.8 x 10 -12 
28 ± 20 x 10 -12 

-3.9F0.7 x 10 -12 
-4.3 ± 1.7 x 10 -12 
37± 19 x 10 -12 

—3.9 ± 0.8 x 10 -12 
-4.4 ± 1.7 x 10 -12 
27± 15 x 10 -12 

-3.9± 1.0 x 10 -12 
-4.1 ± 1.9 x 10 -12 
25± 15 x 10 -12 

-3.5 ±0.8 x 10 -12 
-4.0± 1.7 x 10 -12 
30± 18 x 10 -12 


were assumed to be externally mixed in that the concentra- 
tion of one did not affect the RJF of another. 

The basic quantities we take from these studies are the 
forcing efficiencies (RJF change per unit emission rate change 
) and (for the ozone precursors only) the changes in 

Ac ,™ 4 

methane lifetime per unit emission rate ( — -f - — ). /•", are the 
annual average changes in RF (Wm -2 ) following a change 
in emission rate of species i by E, (kgyr -1 ), and Aa ™ 4 are 
the changes in the methane lifetime (yr). In Sect. 3 the global 
annual average forcing efficiencies from the tables in Yu et 
al. (2013) and Fry et al. (2012) are used; in Sect. 4 we ex- 
pand on the results provided in these papers by splitting the 
forcing efficiencies into 4 latitude bands. 


3 Global metrics 


3.1 Global warming potentials 


GWPs have been the most common way to compare the 
relative climate effects of different species. The absolute 
global warming potential (AGWP) is defined (IPCC, 1990) 
as the integrated RF over a time horizon H due to a 1 kg 
pulse emission of a species and has units Wm - 2 yrkg -1 . 
For a species with an exponential decay with lifetime a, 
AGWP (H) = Aa (l — exp (— H / a )) where A is the specific 
RF (RF per kg atmospheric burden). AGWP is typically nor- 
malised by the AGWP for CO 2 to give a dimensionless quan- 
tity GWP. Here we use the CO 2 impulse response function 
from the IPCC AR4 (Forster et al., 2007) where the RF from 

a pulse of CO 2 is given by A c ° 2 ^ciq ± J2 a; exp ^ ■ 

A c° 2 

is the radiative efficiency for CO 2 in Wm 2 kg 1 
(1.82 x 10 -1 -’’ Wm -2 kg -1 ), a 0-3 are unit less coefficients 
[0.217, 0.259, 0.338 and 0.186], and oq' 0 ^ are CO 2 response 
timescales [172.9, 18.51, 1.186] years. This gives AGWPs 
for CO 2 of 2.5 x 10 -14 and 8.7 x 10 -14 Wm -2 yrkg -1 at 
time horizons of 20 and 100 yr respectively. 


a step change in emission. The equilibrium RFs, F per emis- 
sion change E for each species {FIE) (e.g. Table 6 of Yu 
et al. (2013)), are equal to Aa. Hence, provided that the time 
horizon chosen is much longer than the lifetime of the aerosol 
{a <-£ H), the AGWP is also simply equal to the equilibrium 
response to a unit step change in emissions and is indepen- 
dent of the time horizon. A similar argument was made by 
Bond et al. (2011) where they introduced the specific forcing 
pulse (SFP) (units Jkg -1 ). For short-lived species (« <± //) 
the global SFP is simply a scaling of the AGWP. We do not 
use the SFP concept further in this paper, but it is useful to 
illustrate that for short-lived species the AGWP can be con- 
sidered an effectively instantaneous energy perturbation to 
the system. 

The AGWP is a measure of the globally-averaged net ra- 
diative effect of an emission pulse, but the emissions them- 
selves do not need to be globally emitted. Emissions from 
different regions have different lifetimes and the specific RFs 
may be regionally dependent. Both of these lead to a depen- 
dence of AGWP on the emitting region. AGWPs for SO 2 , 
POM and BC are listed for the 4 HTAP regions in Table 1 ; 
these are the same values as listed in Yu et al. (2013) Ta- 
ble 6 with a change of units. AGWPs for SO 2 and POM 
are negative reflecting the cooling impact. The values for 
EA are slightly smaller than for the other continents due to 
more rapid removal of the aerosols and lower oxidation effi- 
ciency (Koch et al., 2007). Pollutants are typically exported 
from EA by frontal systems with associated strong precipita- 
tion (HTAP, 2010). In contrast, BC emissions have positive 
AGWPs, and the values for EU are slightly larger than for 
the other continents as these absorbing aerosols are advected 
over the brighter surfaces of the Arctic. The GWP values can 
be derived by normalising by the AGWP C ° 2 (see above), giv- 
ing GWP(100)s of -40 ±9, -46 ±20, and 340 ±210 for 
SO 2 , POM and BC, respectively for emissions from “All” 
regions. These are of slightly lower magnitude than in Fu- 
glestvedt et al. ( 2010 ), but within the uncertainty ranges. 

3.1.2 Methane 


3.1.1 Aerosols 

For the aerosol species, the lifetimes are sufficiently short 
that the atmospheric burdens rapidly reach equilibrium after 


The decay time for a methane perturbation is around 12 yr, 
hence for time horizons of 20 or 30 yr, the criterion (a <5jC //) 
no longer holds. For methane, the HTAP experiments were 
based on a global change in the methane concentration 
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rather than regional emission changes to save running the 
chemistry models to equilibrium. Here the direct AGWP is 

F ™4 

again given by — exp(— H / /a™ 4 )) where F™ 4 is 

the methane forcing diagnosed from the methane change 
experiment using the formula from Table 3 of Myhre et 
al. (1998). The factor / accounts for the effect of methane 
on its own lifetime (Prather et al., 1996). The implied 
emissions necessary to cause this change in concentra- 
tion are given by F™ 4 = A 5 CH4 //a CH4 where AF CH4 is 
the change in methane burden. In the HTAP methane per- 
turbation experiment (SR2) where methane was decreased 
by 20% from 1760 ppb to 1408 ppb, AB CH4 =979Tg, 
o! CH4 = 8.7±1.3yr, / = 1.34 ±0.05 (Fiore et al., 2009) giv- 
ing an implied emission of — 84Tgyr _1 . Methane is an 
ozone precursor, it affects the oxidation rate of SO 2 to sul- 
phate aerosol, and it generates water vapour (only impor- 
tant in the stratosphere). The two indirect forcings due to 
changes in ozone (20.5 %± 1.5% of F™ 4 ) and sul- 
phate F<?^ 4 (— 1 . 1 % ± 2.4 % of Fch 4 4 ) were calculated in Fry 
et al. (2012); we also assume a stratospheric water vapour 
forcing ec l ua l t0 15% ±10% of the direct methane 

forcing (Forster et al., 2007) giving a total indirect contribu- 
tion of 34 ± 10 %. Contributions from ozone changes above 
the tropopause are not include here (see Sect. 5.1). The total 
direct ± indirect AGWP from methane (excluding CO 2 ) is 
therefore 

AGWP™ 4 = (1) 

/ Z 7 CH 4 , £-03 , ]r.S04 , r H20, £i / [j / 4-„.CH 4 \, 

™ch 4 + F CH 4 + F ch 4 + F CH 4 > A gCH 4 (1 - exp [~H/ fa 4 j) 

3.1.3 Ozone precursors 

Ozone precursors (CO, VOCs and NO x ) have a short-term 
impact on RF by generating ozone and perturbing sulphate 
oxidation; they also have a longer-term impact by chang- 
ing the methane lifetime. The short-term impacts on ozone 
F ( ° 3 and sulphate F ( S ° 4 forcing are diagnosed from the equi- 
librium experiments perturbing emissions of species i but 
holding methane fixed. For H a, as for the aerosols, the 
AGWPs for the short-lived components are then simply 
(F)° 3 ± F- S ° 4 )/F/. Although methane concentrations were 
fixed in these experiments, the longer-term impacts on 
methane can be diagnosed from the change in the fluxes 
through the CH 4 ± OH reaction in the models. Changes in 
methane lifetime A^ 114 are then calculated from the flux di- 
agnostics. If the methane concentrations had not been fixed 
the methane burden would have responded to the changes 
in lifetime moving towards a new equilibrium burden of 

( A " C > C " 4 ) B™ 4 where F CH4 is the unperturbed burden 
(4983 Tg in the HTAP study). For small perturbations, the 
impact of these methane lifetime changes on the methane 
burden can be shown to be equivalent to the impact of 


1 • • An -™ 4 R ™ 4 o- 1 

methane emissions ot “ CH4 Since the response to 

methane emissions has already been derived (AGWP (I ' 4 ), 
the AGWP for ozone precursors including the short and long- 
term (methane related) components is given by 

AGWP 1 ' = -^r (f° 3 ± F S ° 4 j ± x , AGWP™ 4 (//), (2) 

where the scaling factor Xi = represents the 

equivalent methane emission and the AGWP™ 4 implicitly 
includes the long-term ozone and sulphate responses. Note 
only the long-term components are functions of the time 
horizon. The GWPs at 20 and 100 yr for the ozone precur- 
sors were presented by Fry et al. (2012), their Table S2. 


3.2 Global temperature-change potentials 


The GTP metric (Shine et al., 2005) is more closely re- 
lated to observable climate impacts than the GWP, as it is 
calculated from temperature changes. The absolute global 
temperature-change potential (AGTP) is given by the surface 
temperature response following an instantaneous pulse emis- 
sion. The global surface temperature responds to changes in 
the RF at the tropopause on a spectrum of timescales. We 
follow the same methodology in Fuglestvedt et al. (2010) 
and use a temperature response function that is the sum of 

two exponentials R(t) = % ex P where the coeffi- 

cients cj and dj are the climate sensitivities and timescales 
of the two modes. From Boucher and Reddy (2008) 
we use ci =0.631 K(Wm- 2 )- ! ), c 2 = 0.429 K(Wm- 2 )-*), 
di = 8.4 yr, d 2 = 409.5 yr, giving an equilibrium climate sensi- 
tivity of 1.060 K(Wm -2 ) -1 ) (3.9 K for C0 2 doubling). The 
GTPs are constructed by normalising by the corresponding 
AGTP for CO 2 . For this we use the expression in Fuglestvedt 
et al. (2010) (their Eq. A3). 


AGTP™ 2 (// ) = A™ 2 


j= 1 


aocj I 1 — exp 




± ( 3 ) 


Using the CO 2 radiative efficiency (A c ° 2 ) and the response 
parameters a,- and af° 2 as described in Sect. 3.1 gives 
AGTP™ 2 of 6.8 x 10 -16 and 5.1 x lO-^Rkg" 1 for 20 and 
100 yr time horizons respectively. 


3.2.1 Aerosols 


For aerosols with lifetimes much less than the time horizons 
of interest and the shortest climate timescale (a <^H, a <£d\) 
the AGTP is simply a scaling of the temperature response 
function by the equilibrium forcing from a unit step change 
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Table 2. GTPs at 20 and lOOyr time horizons for the HTAP aerosols (or precursors), methane and non-methane ozone precursors. Data are 
given for emissions from the 4 HTAP continents separately and combined. For methane only a global number is given. Uncertainties are inter- 
model differences in the CTM responses to emission perturbations from Yu et al. (2013) and Fry et al. (2012). This includes uncertainties in 
the aerosol and ozone responses, and in the methane lifetime response. Note these values exclude the effects of aerosol cloud interactions, 
the effect of BC deposition on snow albedo, changes in ozone above the tropopause and carbon-cycle interactions as described in Sect. 5. 



EA 

EU 

NA 

SA 

All 

GTP(20) 

so 2 

-31 ± 13 

—43 ±15 

-43 ± 14 

-42 ±9 

—38 ± 12 

POM 

-50± 17 

—58 ±18 

—59 ±15 

—58 ±17 

—55 ± 16 

BC 

410 ±350 

530 ± 190 

400 ± 150 

350 ± 160 

420 ± 190 

ch 4 





54.6 ±8.8 

NO x 

—55.6 ± 23.8 

-48.0 ± 14.9 

—61.9 ± 27.8 

-124.6 ±67.4 

-62.1 ±26.2 

NMVOC 

8.4 ±4.6 

9.5 ±6.5 

8.6 ±6.4 

15.7 ±5.0 

10.0 ±5.7 

CO 

3.5 ±1.3 

3.2 ±1.2 

3.7 ±1.3 

3.4 ± 1.0 

3.5 ±1.2 

GTP(100) 

so 2 

-4.1 ± 1.6 

—5.7 ±2.0 

—5.7 ±1.8 

—5.5 ±1.1 

-5.1 ± 1.6 

POM 

-6.7 ±2.0 

-7.7 ±2.4 

-7.8 ±2.0 

-7.7 ±2.3 

—7.3 ±2.1 

BC 

55 ±46 

71 ±25 

53 ±20 

46 ± 21 

56 ±25 

ch 4 





3.6 ±1.2 

NO x 

— 1.3 ± 2.1 

—2.5 ±1.3 

— 1.7 ±2.1 

—4.6 ±5.1 

—2.2 ±2.1 

NMVOC 

0.7 ±0.4 

0.8 ±0.5 

0.7 ±0.5 

1.3 ± 0.5 

0.9 ±0.5 

CO 

0.26 ± 0.12 

0.24 ±0.11 

0.27 ± 0.12 

0.27 ± 0.10 

0.26 ±0.11 


in emission. 

AGTP' (H) = —R(H) = AGWP' x R(H) 

Ei 

The GTP is calculated by dividing this by AGTP C ° 2 . The 
GTPs for the 3 HTAP aerosol species from each of the 4 
continents and combined (“All”) are shown in Fig. 1 for 
a 20 yr time horizon. The values for a lOOyr time hori- 
zon would have exactly the same pattern, but scaled by 

t^ioo) AGTP ~ (20> (=0.138). GTPs for both time horizons 
2? (20) AGTP co 2(100) v ' 

are listed in Table 2. The patterns are also the same as for the 
AGWPs in Table 1 and the forcing efficiencies in Table 6 of 
Yuet al. (2013). 

3.2.2 Methane 

For methane, the interactions between the methane lifetime 
and the climate response timescales become important. Com- 
bining the expression for the methane AGWP with equation 
A4 from Fuglestvedt et al. (2010) we get for the total AGTP: 

AGTP™ 4 (H) = (C™ 4 + F c % + F C S ° 4 + f"£) • (4) 

g7i^(“ p (7i^)- exp (w)) 

The GTP™ 4 for 20 and lOOyr are shown in Fig. 2 and are 
listed in Table 2. The decrease in GTP between 20 and 1 00 yr 
is a factor of 15, much larger than for the aerosols, This is be- 
cause the AGTP™ 4 has not dropped off as quickly by year 


20 as it did for the aerosols. The dominant contribution to the 
QjpCH 4 Llncer t a inty comes from the inter-model variability 
in the methane lifetime in the control run (a™ 4 ). The di- 
rect methane contribution to the GTP(20) and GTP(100) are 
in agreement with Boucher et al. (2009), however the ozone 
contribution is slightly lower here as discussed in Sect. 5.1. 

3.2.3 Ozone precursors 

The AGTPs for the other ozone precursor species can be cal- 
culated by combining the equilibrium responses with the ef- 
fects on the methane lifetime to obtain an analogous expres- 
sion to the AGWP: 

AGTP'' (H) = — (f° 2 + Ff° 4 ) R(H) + X i AGTP™ 4 (H) (5) 

The time evolution of the total AGTPs therefore depends on 
the relative balance between the short and long-term compo- 
nents, i.e. the almost instantaneous change in ozone and the 
decadal timescale for changes in methane. Figure 3 shows 
the GTP for NO x , NMVOC and CO emissions (summed over 
the 4 continents). For NO x emissions the long-term compo- 
nent acts in the opposite sense to the short-term because NO x 
emissions produce ozone, but decrease methane (/, is neg- 
ative). The AGTP starts positive, but becomes negative af- 
ter 10 yr. For NMVOC and CO emissions the short-term and 
long-term components have the same sign (as they increase 
both ozone and methane), with the long-term component be- 
ing relatively more important for CO. In all cases, the GTP 
is dominated by the long-term component from about 30 yr 
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so 2 



EA EU NA SA All 

Emission region 



EA EU NA SA All EA EU NA SA All 

Emission region Emission region 


Fig. 1. GTP(20) for emissions of SO 2 , POM and BC from the four HTAP continents. Error bars reflect 1 s.d. uncertainty in the multi-model 
responses. 


Table 3. Regional response coefficients in K(Wm -2 ) -1 for unit forcing from changes in BC and O 3 . Values are taken from Fig. 1 of Shindell 
and Faluvegi (2009). Data from CO 2 changes are used for the 90° S-28° S forcing region. The k] m coefficients used in section 4 are derived 
by normalising these values by the climate sensitivity. 


BC 





Forcing region 




90° 

S-28° S 

28° 

S-28° N 

28° N-60° N 

O 

O 

50 

N-90° N 



90° S-28° S 


0.19 


0.06 

0.02 


0.00 

Response region 

28° S-28° N 


0.09 


0.17 

0.07 


0.02 

28° N-60° N 


0.07 


0.24 

0.14 


0.08 


60° N-90° N 


0.06 


0.31 

0.15 


-0.08 

O 3 





Forcing region 




90° 

S-28° S 

28° 

S-28° N 

28° N-60° N 

O 

O 

50 

N-90° N 


90° S-28° S 


0.19 


0.13 

-0.06 


-0.03 

Response region 

28° S-28° N 


0.09 


0.26 

0.09 


0.02 

28° N-60° N 


0.07 


0.15 

0.20 


0.06 


60° N-90° N 


0.06 


0.13 

0.05 


0.07 


to 60 yr. Beyond 60 yr the short and long-term components 
are of similar magnitudes, although both are very small. The 
results of this can be seen in the bar charts of Fig. 4 and val- 
ues listed in Table 2 which break the GTPs down by emission 
species and emission region. The net NO x GTPs are negative 
with the largest impact coming from emissions in SA, and 
the least impact from EU. For VOCs the largest impact is also 
from SA. The tropics are the region of greatest methane oxi- 
dation, so NO x and VOC emissions from SA have the largest 
impact on methane. CO has a sufficiently long lifetime that 
it can be advected away from the source continent; hence, 
the climate impact is largely independent of emission region. 
The fractional contribution to the GTP from the methane re- 
sponse is greater for the 20 yr time horizon than the lOOyr 
horizon. Although counterintuitive, the short-term contribu- 
tions have already fallen rapidly by the 20 yr horizon and 
the curves (Fig. 3) begin to flatten (decreases of a factor of 
7 between 20 and lOOyr), whereas for the long-term con- 
tributions (driven by methane changes) the 20 yr horizon is 
roughly in the steepest portion of the curve, declining there- 
after by a factor of 15 by the 100 yr horizon. 

The error bars are a combination of the uncertainties 
in /f° 3 , F S ° 4 , a™ 4 , and Act™ 4 from the inter-model 
variability added in quadrature. The dominant contribution 


comes from the inter-model variability in the changes to the 
methane lifetime (Aa™ 4 ). 

4 Regional climate change 

Shindell and Faluvegi (2009), Shindell and Faluvegi (2010) 
and Shindell (2012) introduced the idea of the absolute re- 
gional temperature-change potential (ARTP) to quantify re- 
gional climate change from heterogeneous forcing patterns. 
There is not a one-to-one correspondence between the forc- 
ing in a region and the local temperature response, as there 
are RF contributions from outside the regions due to the 
transport of heat. Shindell and Faluvegi (2009) showed that 
the forcing-temperature relationship could be characterised 
using 4 broad latitude bands: southern mid-high latitudes 
(90° S-28° S), tropics (28° S-28° N), northern mid-latitudes 
(28° N-60° N), and the Arctic (60° N-90° N). Thus the ex- 
pression for the AGTP can be generalised as: 

H 

ARTP m (H) = L/ ki m %P-R(H -t)dt where ARTP m is 
l 0 

the temperature response in latitude band m, lq m are the re- 
sponse coefficients derived in Fig. 1 of Shindell and Faluvegi 
(2009) normalised by the climate sensitivity, and F[ is the 
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ch 4 



Global 


Fig. 2. GTP(20) and GTP(IOO) for global emissions of CH 4 . The 
GTPs are split according to the contributing components. The con- 
tribution from stratospheric water vapour is included in with CH 4 . 
Error bars reflect 1 standard deviation uncertainty in the multi- 
model responses. 


tropopause RF in latitude band / for a pulse emission E. For 
SO4, POM and CF4 we use the average k\ m for CO2 and 
SO4 (Shindell and Faluvegi, 2010). For BC and O3 we use 
the ki, n appropriate for the species except for responses to 
forcing in the 28° S-90° S band which was only derived for 
CO2 by Shindell and Faluvegi (2009). The coefficients used 
for BC and O3 are taken from Table 3, normalising by the 
climate sensitivity. We assume that the temperature response 
function R applies equally to temperature changes in each 
region. As discussed in Shindell and Faluvegi (2010), in the 
different latitude bands encompassing such diverse regions 
as the Arctic and the tropics, the partitioning between the 
fast and slow responses and even the time constants of these 
responses may vary from the global average. There is no data 
currently available to break this down further, although Shin- 
dell (2012) suggests that the land response might be 20% 
larger than the average of land and ocean. 

4.1 Aerosols 

For aerosol emissions, the ARTPs have the same functional 
form as the AGTPs with the magnitudes given by the prod- 
uct of the response coefficients with the forcing pattern. The 
temperature change in latitude band m due to the emission of 


species i : 

ARTP;„(7/) = V ^-kimAGTP 1 (H) 

I Fi 

where k) is the equilibrium globally-average forcing as be- 
fore, and Fij is the RF averaged over latitude band /. The lat- 
itudinal patterns of the ARTPs for aerosols do not vary with 
time. Figure 5 shows the normalised equilibrium forcing pat- 
terns ( Fjj/Ej ) and the ARTPJ n (20) values for the 3 aerosol 
species. Note that for the scattering aerosols (SO2 and POM 
emissions) the largest RF is in the latitude band closest to the 
emitting region (28° N-60° N for EA, EU and NA emissions, 
28° S-28° N for SA emissions). There seems to be very lit- 
tle RF from scattering aerosols over the bright Arctic lati- 
tudes as there is less contrast between the brightness of the 
aerosol and the underlying albedo. For the absorbing aerosol 
(BC), due to the high contrast, the RF is high over the Arc- 
tic latitudes, and for EU emissions the Arctic RF exceeds 
that of the northern mid- latitudes. Applying the matrix of re- 
sponse coefficients smoothes out the impacts such that there 
is less latitudinal variation in the temperature responses than 
in the RFs. The Arctic temperatures are particularly sensitive 
to forcing changes in the northern mid-latitudes (at least in 
the coefficient set used here), where the Arctic responses to 
the scattering aerosols are comparable to the northern mid- 
latitude responses. 

For BC the Arctic shows a lower surface temperature re- 
sponse than would be expected from the RF pattern. This is 
because for BC the ki m are negative for the Arctic response 
to Arctic forcing (see Table 3), i.e. a positive RF over the 
Arctic due to changes in atmospheric BC has a cooling ef- 
fect on Arctic temperatures (figure Id of Shindell and Falu- 
vegi 2009). This relationship, while surprising, has been con- 
firmed by other studies (Sand et al., 2013; Flanner, 2013) and 
is caused primarily by a decrease in the meridional heat trans- 
port. The Arctic cooling effect only applies to BC changes 
in the mid to upper Arctic troposphere, such as found here 
following advection from northern mid-latitude and tropical 
sources. Local Arctic emission of BC has a warming effect, 
as does deposition of BC on to Arctic snow from any source 
region (Flanner, 2013). The net effect on the Arctic of BC 
emissions is still warming from each of the HTAP regions 
since RF changes in the northern mid-latitudes and tropics 
have a large positive impact on Arctic temperatures by in- 
creasing the meridional heat transport. In fact in this study 
BC emissions from SA have the largest impact on Arctic 
temperatures since they exert a strong RF over the tropics 
and northern mid-latitudes and little over the Arctic. The re- 
sponses in the Arctic to BC emissions might be expected to 
be more sensitive to the more northerly source regions ifBC 
deposition to snow were accounted for. 

The ARTP(20) values are shown, but the latitudinal pat- 
terns will be identical for any time horizon, albeit with dif- 
ferent overall scaling. The aerosol impacts on the northern 
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Fig. 3. GTPs as functions of time horizon for emissions of NO x , NMVOC and CO from all 4 continents combined. The short and long-term 
components are also shown. 


NO, 


voc 


CO 




Fig. 4. GTP(20) and GTP(IOO) for emissions of NO x , VOC and CO from the four HTAP continents separately and combined (“All”). The 
GTPs are split according to the contributing components. For O 3 and SO 4 these combine both the short and long-term components. The 
contribution from stratospheric water vapour is included in with CH 4 . Error bars reflect 1 standard deviation uncertainty in the multi -model 
responses. 


latitudes are up to twice that implied just using the global av- 
erage response. The “global” ARTP bars are calculated from 
an area-weighted average of the responses in the 4 latitude 
bands. These are slightly larger than the calculated AGTPs in 
the previous section, which implicitly used global averaged 
climate sensitivity. This difference can be considered as an 
efficacy larger than unity due to the forcing being concen- 
trated at the higher northern latitudes where the climate re- 
sponse is stronger. 


4.2 Methane 

Analogous to the aerosols, for methane 
ARTP ™ 4 (H) = kim AGTP ™ 4 {H) 

where F C \ U j = F™ 4 , + Fq Ul + Fquj’ notin S that the di- 
rect methane forcing and indirect forcings through ozone and 
sulphate will have different latitudinal distributions. We do 
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Fig. 5. Normalised equilibrium forcing in mW m 2 Tg 1 (upper panel) and ARTP(20) in mKTg 1 (lower panel) for emissions of SO 2 , 
POM and BC from the four HTAP continents. The coloured columns represent the forcing and temperature response in 4 different latitude 
bands (90° S-28° S, 28° S-28° N, 28° N-60° N, 60° N-90° N) and globally. 


not include the stratospheric water vapour forcing here as 
we have no information on its latitudinal profile. The RF 
and ARTP distributions are shown in Fig. 6. For methane, 
the forcings are more evenly distributed as methane is well- 
mixed, with the maximum RF in the tropics, but a signif- 
icant RF in both the southern and northern extra-tropics. 
The ARTPs show less temperature response in the southern 
extra-tropics as this region is less sensitive to RF changes 
outside these latitudes. As for the aerosols, the latitudinal pat- 
tern of the temperature response to methane emissions does 
not vary with time. 

4.3 Ozone precursors 



For non-methane ozone precursors, the ARTPs are a combi- 
nation of two functional forms: 


ARTP;„(H) = 

E (*/« (^°/ 3 + F u 4 )) R W + x/ARTp™4 ( //) 

Here because the latitudinal patterns of the Fj j and Fch 4 .i are 
different the overall spatial pattern of the ARTP does evolve 
over time. 

Figure 7 shows the latitudinal components of the forc- 
ing patterns for each ozone precursor from each region. 
The top panel shows the normalised short-lived equilibrium 
forcing terms (F° 3 + F^ 4 ) j The lower panel shows 
the total equilibrium forcing patterns including the long- 
term component due to changes in the methane lifetime 
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Fig. 6 . Normalised equilibrium forcing in mW nr Tg 1 (left), and 
ARTPs 20 and 100 (right), for global emissions of CH 4 . The 
coloured columns represent the forcing and temperature responses 
in 4 different latitude bands (90° S-28° S, 28° S-28° N, 28° N- 
60° N, 60° N-90° N) and globally. 
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Fig. 7. Short-term (upper) and total (lower) normalised equilibrium forcing in mW m 2 Tg 1 for emissions of NO x , VOC and CO from the 
four HTAP continents. The coloured columns represent the forcing in 4 different latitude bands (90° S-28° S, 28° S-28° N, 28° N-60° N, 
60° N-90° N) and globally. 


(7}°/ + F^ 4 ) j Ej + xi Fch 4 ,/ along with the forcing pattern 
for methane itself. The RFs depend on the latitude of the 
emissions, with the highest short-term forcing in the north- 
ern mid-latitudes for emissions from EA, EU and NA, and 
in the tropics for emissions from SA. For NO x , the /, are 
negative (reduced methane) so the net forcing is the differ- 
ence between short and long-term components. The net forc- 
ings are negative in the Southern Hemisphere and tropics 
where methane destruction outweighs the ozone production, 
but positive in the northern mid and high latitudes (except 
for SA emissions which have negative forcing in the north- 
ern high latitudes). For NMVOC and CO emissions, the Xi 
are positive (increased methane) so the long and short-term 
components add. When adding the smoother methane terms, 
the latitudinal distributions for the total forcing become more 
even than for the short-term components. For CO, the highest 
forcing is in the tropics for all emission regions. 

The corresponding ARTPs are shown in Fig. 8. For the 
ARTP(20)s, the long-term component is more important 
than the short-term component and the ARTPs for the non- 
methane ozone precursors look similar to the methane ARTP, 
with an opposite sign for NO x emissions. The highest re- 
sponse is in the northern mid-latitudes for VOCs and CO. For 
the ARTP(IOO) the contributions of the short and long-term 
components are more balanced. ForNO x the ARTP(100)s are 
close to zero in the northern mid-latitudes for emissions from 
EA and NA reflecting the warming from ozone balancing the 
cooling from decreased methane. 

All the results from HTAP were based on changing North- 
ern Hemisphere emissions. We do not have results from 


Southern Hemisphere experiments, but might expect the 
forcing and temperature responses to be stronger in the hemi- 
sphere in to which the species are emitted. 

5 Limitations and uncertainties 

We acknowledge limitations to this study in this section, 
most importantly the neglect of the aerosol-cloud interac- 
tions. As a consequence, care must be taken in using some 
of the numerical values presented in Table 2. 

5.1 Changes in stratospheric ozone 

As previously discussed in Fry et al. (2012) the ozone forc- 
ing here only includes ozone changes below the tropopause. 
Sovde et al. (2012) calculate a forcing contribution from 
ozone responses above the tropopause of 15% of the total. 
Adding this would take the ozone contribution to the methane 
metrics to 24 % of the direct methane effect. This ends up be- 
ing a similar value to the 25 ± 15 % often used (e.g. Forster et 
al., 2007) which came originally from the troposphere-only 
responses of a mixture of 2-D and 3-D models discussed in 
chapter 2 of IPCC (1994). 

5.2 Aerosol interactions with clouds and snow albedo 

The responses to aerosol emission changes were calculated 
using chemical transport models that omit interactions with 
clouds, so the RFs presented here only include direct effects 
based on the aerosol column loadings (Yu et al., 2013). The 
indirect effects of aerosols are likely to be as great as the 
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Fig. 8. ARTP(20) and ARTP(100) in mK/Tg for emissions of NO x , VOC and CO from the four HTAP continents. The coloured columns 
represent the temperature responses in 4 different latitude bands (90° S-28° S, 28° S-28° N, 28° N-60° N, 60° N-90° N) and globally. 


direct effect for hydrophilic aerosols (sulphate and POM in 
this study) and may also be important for BC. The indirect 
effects are likely to vary with emission region depending on 
their proximity to stratiform clouds and the regional aerosol 
loading. The impact ofBC deposition on snow albedo is also 
very uncertain. Bond et al. (201 1) estimate that it adds about 
20% to the RF for BC emissions from NA, EU and SA, 
and about 10 % for EA emissions. This might be expected to 
preferentially increase the temperature responses in the Arc- 
tic latitudes. 

5.3 Nitrate aerosols 

NO x emissions have an additional cooling effect by generat- 
ing ammonium nitrate aerosols in regions of high ammonia 
abundance. Unfortunately, nitrate aerosols were not included 
in many of the HTAP simulations, but we can use the results 
of Bauer et al. (2007) who calculated a nonnalised direct RF 
from global anthropogenic NO x emissions of -2.0 x 10“ 12 
Wm _2 kg _1 . This would add —25 and —3.3 to the “All” NO x 
GTP(20) and (100) and (—62 and —2.2 in Table 2), making 
these consistently negative. 

5.4 CO 2 from methane oxidation 

Methane is oxidised in the atmosphere, with a molar frac- 
tion s of between 0.5 and 1.0 resulting in CO 2 (Boucher 
et al., 2009). This provides an emission of CO 2 equal to 
£ x (44/16) x [CH^/a™ 4 which has its own climate impact 
(the factor 44/16 converts from a molar fraction to a mass 


fraction). Thus the addition to the methane GTP is: 

AAGTP^ ation = f ex^x 

0 

^55 “ P (" Tisc) AGTpC °’ (H-I)dl 

The methane oxidation adds 1.3 ±0.4 and 2.1 ±0.7 to 
the GTP(20) and GTP(100) respectively as in Boucher et 
al. (2009). Note that for non-fossil methane this extra CO 2 
is still produced as described above, but should be compen- 
sated for by previous uptake of CO 2 . Thus the above cor- 
rection should not be applied to non-fossil methane without 
also accounting for the uptake, requiring detailed life-cycle 
analysis (for instance of the production of cattle fodder). 

5.5 Carbon cycle effects 

Changes to the atmospheric composition of gases and 
aerosols can also affect the carbon cycle thus adding further 
indirect climate impacts (e.g. Sitch et al., 2007). Collins et 
al. (2010) showed that adding in the effects of surface ozone 
on vegetation damage and reduced uptake of CO 2 added 
about 10% to the methane GTP(20) and could change the 
sign of the NO x GTPs. Hence the sign of the net impact of 
NO x emissions is still uncertain. 

The temperature changes calculated in the GTPs can also 
lead to an impact on the carbon cycle, for instance through 
increased soil respiration. The above calculations of GWP 
and GTP for species i as ratios AGWP'/AGWP C ° 2 and 
AGTP'/AGTP C ° 2 treat this inconsistently (as do many other 
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Fig. 9. The direct methane GTP (excluding ozone, sulphate and 
stratospheric water vapour), and additional contributions from CO2 
from two sources: the atmospheric oxidation of methane; and the 
impact of the methane-induced temperature change on the carbon 
cycle. 

studies, e.g. Forster et al., 2007). The temperature-carbon 
cycle feedback is implicitly included in the response func- 
tion used for CO 2 (Joos et al., 2001) and hence in the 
denominator AGTP C ° 2 , but is excluded in the numerator 
AGTP'. Gillett and Matthews (2010) and Sarofim (2012) 
using carbon cycle models find that this can underestimate 
the GTPs for methane by 20 %. A way of obtaining a sim- 
ple estimate of this effect is to make the crude assumption 
that the change in the land-atmosphere flux of CO 2 is lin- 
early proportionally to the surface temperature change with 
a coefficient T. Multi-model comparisons of carbon cycle 
models (Friedlingstein et al., 2006; Arora et al., 2013) sug- 
gest a change in the land and ocean uptake of approximately 
— 1 GtCyr -1 per K temperature change. We therefore use 
a T of +(44/12) x 10 12 kg(CC>2)yr _1 K -1 . We associate an 
uncertainty of 100% with this to reflect a range from no ef- 
fect to twice what we have quoted. Treating the decreased 
uptake as an extra CO 2 release of F x AGTP' gives an addi- 
tion to the AGTP due to the climate-carbon cycle feedback 
of: 

H 

AAGTP' Hmate = / r x AGTP' (f) x AGTP C ° 2 ( H - t)dt 
0 

Note that because the AGTP C ° 2 does include the climate 
feedback we do not need to iterate this calculation any fur- 
ther. The contributions to the methane GTP from methane 
oxidation and climate-carbon cycle feedback are illustrated 
in Fig. 9 along with the direct methane impact (i.e. exclud- 
ing ozone, sulphate and stratospheric water vapour). The 


climate-carbon cycle feedback effects are potentially large 
adding 1 .7 ± 1 .7 and 4.3 ± 4.3 to the methane GTPs at 20 and 
100 yr (55 and 3.6 in Table 2). Note that the climate-carbon 
cycle feedback will also affect GWPs too, adding 4.5 ±4.5 
to the methane GWP(100) value (24.2 in Fry et al., 2012). 
For aerosols and the short-term component of ozone precur- 
sors this climate-carbon cycle feedback could add 17% and 
170% to the 20 and 100 yr GTPs and 4% and 20% to the 
20 and 100 yr GWPs. The CO 2 impacts are more important 
over longer timescales as they persist for considerably longer 
than the direct effects, since CO 2 continues to be released for 
some time after the pulse of the original emitted species has 
decayed. 

5.6 Temperature response 

One major uncertainty is in the assumptions about the tem- 
perature response, both the timescales of the temperature re- 
sponse function R and the latitudinal patterns included in the 
k /„,. Oliie and Peters (2012) explore the sensitivity of GTP 
to model spread in R and find that uncertainties increase 
with time horizon and decrease with lifetime. For GTP™ 4 an 
uncertainty of — 12/+10% was found for a 20 yr time hori- 
zon, which increased to —11 1+15 % for a time horizon of 
100 yr (for 5th— 95th percentiles). Fligher uncertainty ranges 
were calculated for short-lived species such as aerosols; for 
20 yr time horizon these were — 59/+85% and — 88/-+80 % 
for 100 yr. 


6 Conclusions 

We have examined the climate impacts of aerosols and reac- 
tive gases, characterising them using metrics to inform policy 
makers of the relative importance of the different species, the 
dependence of impacts on emissions location, and the poten- 
tial impacts of control policies. The analytical formulae con- 
structed here could be used to calculate metrics for similar 
experiments in the literature (or yet to be performed) that re- 
port RFs and methane lifetime changes. 

For the aerosol species the chain from RF to GWP to GTP 
adds little extra information on the relative differences since 
they are simply scalings of each other. Hence the relative 
global importance of emissions of different species from dif- 
ferent regions remains unchanged from Yu et al. (2013). The 
inter-model variability is larger than the variability between 
emission regions, suggesting that a single metric averaged 
over all emitting regions would be a usable approximation. 
The small differences between regions that are suggested are 
a lower climate impact for the emission of SO 2 and POM 
from East Asia, and a higher climate impact for the emission 
of BC from Europe. 

We calculate GTP CH4 s of 56 ±9 and 3.7 ±1.2 for 20 
and lOOyr after including a correction for the expected in- 
creases in lower stratospheric ozone. Adding in the effects of 
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climate-carbon cycle feedbacks and CO2 from methane oxi- 
dation would increase these to 59 ± 9 and 10. 1 ± 4.5. 

For ozone precursors, the relative importance of the short 
and long-term components varies with the time horizon used. 
This is most noticeable for NO x emissions where the short 
and long-term components act in opposite directions so the 
GTpNO*s are dominated by the negative methane component 
on the 20 yr timescale, but after 100 yr the ozone compo- 
nent plays a more significant role. As with the GWPs (Fry 
et al., 2012), the GTPs for the ozone precursors are largely 
independent of source region for emissions in the north- 
ern mid-latitudes, with larger magnitudes for NO x and VOC 
emissions from the more tropical SA region due to more ac- 
tive photochemistry. 

The different patterns of responses when moving down the 
impact chain become more apparent when breaking down 
the temperature changes into latitude bands. There is no di- 
rect relationship between the RF in a region and the surface 
temperature change in that region. Even when using broad 
latitude bands, there are contributions to the surface tem- 
perature from RF outside the region, although that contribu- 
tion decreases with latitudinal distance. Thus, the latitudinal 
temperature patterns can differ significantly from the global 
mean, and also from the latitudinal forcing patterns. 

The RF from sulphate and POM aerosols peaks strongly 
over the latitudes of the emissions, and so the northern mid- 
latitude temperature response to northern mid-latitude emis- 
sions is approximately twice as large as the global average. 
For BC there is a large RF component over the Arctic; how- 
ever this leads to a cooling of the Artie surface. Therefore 
the Arctic temperatures are most sensitive to emissions from 
SA since a larger proportion of BC emitted from the other 
more northerly regions is advected to the Arctic upper tropo- 
sphere. This excludes the effects on Artie temperatures ex- 
pected from BC deposition to snow which might increase 
the importance of the more northerly emission regions on the 
Arctic. 

The RFs and RTPs for CFI4, VOC and CO emissions show 
smoother latitudinal profiles than for the aerosols with a typ- 
ically 20-30 % greater temperature response in the northern 
mid-latitudes than the global average. NO x emissions from 
the Northern Hemisphere source regions lead to negative RFs 
in the Southern Hemisphere and the tropics and positive RFs 
in the northern mid-high latitudes. The resulting temperature 
changes however are cooling in all latitudes, with the excep- 
tions of northern mid-latitude responses to EA and NA emis- 
sions on the 100 yr timescale (although the lOOyr responses 
are extremely small). 

The findings in this paper begin to address important pol- 
icy questions of whether emissions controls within a region 
have a larger climate effect in that region, and whether NTCF 
mitigation controls can be optimised for their impacts on vul- 
nerable regions such as the Arctic. 
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